Two-orbital quantum spin model of magnetism in the iron pnictides 
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We study a two-orbital spin model to describe (tt, 0) stripe antiferromagnetism in the iron pnic- 
tides. The "double-spin" model has an on-site Hunds's coupling and inter-site interactions extending 
to second neighbors (inter- and intra-orbital) on the square lattice. Using a variational method based 
on a cluster decomposition, we optimize wave functions with up to 8 cluster sites (up to 2^® vari- 
ational parameters). We focus on the anomalously small ordered moments in the stripe state of 
the pnictides. To account for it, and large variations among different compounds, we show that the 
second-neighbor cross-orbital exchange constant should be ferromagnetic, which leads to "partially 
hidden" stripe order, with a moment that can be varied over a large range by small changes in 
the coupling constants. In a different parameter region, we confirm the existence of a canted state 
previously found in spin- wave theory. We also identify several other phases of the model. 

PACS numbers: 74.70.Xa, 75.10.Jm, 75.50.Ee, 75.30.Et 



I. INTRODUCTION 

Like the cuprate high-Tc superconductors, the more 
recently discovered iron pnictide superconductorsii^ also 
exhibit interesting magnetic propertiesi^"— While the 
antiferromagnetism of the Mott-insulating parent com- 
pounds of the cuprates can be very well accounted for by 
the two-dimensional (2D) S = 1/2 Hcisenbcrg modclfS 
there is no such simple reference point for the Fc pnic- 
tides, in which the electrons are always itinerant. Quan- 
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tum spin models have nevertheless been employecti 
describe some of their intriguing magnetic properties, as 
more manageable alternatives to t-J^ or Hubbar d^^i^^ 
models of itinerant electrons. In contrast to the (tt, tt) 
(Neel) order of the cuprates, most of the pnictides ex- 
hibit (tt, 0) "stripe" order, which can be achieved with 
frustrated spin modelsJ^ An on-going theoretical chal- 
lenge is to explaiufiii^ in a general sense, the ordered 
moments, which are often anomalously small but vary 
considerably among different compounds i^^— 

We here address this issue, as well as other generic as- 
pects of multi-orbital magnetism, using a model with two 
electronic spins per lattice site. While the neglect of the 
charge degrees of freedom will clearly also affect the mag- 
netic properties to some extent, it is still interesting and 
important to ask the question of whether a correct effec- 
tive low-energy description of the magnetic properties of 
the Fe pnictides is possible within such a simplified spin- 
only model. The two-orbital spin system should be much 
more realistic in this regard than the often used models 
with single S* = 1/2 or 5 = 1 spins on the sitesj^ii^ii^ We 
will include in the two-spin description an on-site Hund's 
coupling, as well as inter- and intra-orbital exchange be- 
tween nearest and ncxt-ncarcst-ncighbor spins. 

Studying this highly frustrated quantum spin system 
with a cluster-mean-field method (a variational cluster- 
product-state ansatz where the wave function on small 
clusters is optimized fully), we find a very rich phase di- 
agram, including a regime in which the behavior seems 



appropriate for describing the small-moment stripe anti- 
ferromagnetism of the pnictides. 



II. DOUBLE-SPIN MODEL 

The valence electrons in the Fe pnictides reside pre- 
dominantly in Fe d orbitals. In a simplified picture suit- 
able for some compounds, the low-energy bands^ can be 
reproduced using only d^z and dyz orbitals^i (although 
in some cases other orbitals may also be neede d^ ^'^^ ). 
In undoped systems there are two electrons per Fe in 
this case. Despite the electron itineracy, it has been 
argued that the magnetic properties are still correctly 
maintained upon further reducing the two-orbital Hub- 
bard model to a "double-spin" model (DSM)fiii^ where 
two localized 5=1/2 spin degrees of freedom per lattice 
site are retained. The "A" and "B" spins are coupled to 
each other intra-site by a ferromagnetic Hund's coupling 
Jh, and significant inter-site couplings should extend up 
to second-nearest neighbors on the square lattice of iron 
ions. These interactions originate from supcrexchange 
through As ionsi^I 

The spin-only description of course neglects some im- 
portant aspects of the pnictides. One appealing proposal 
to explain their small and highly varying magnetic mo- 
ments is a "semi-localized" picture, in which the elec- 
trons on one orbital are essentially localized, while those 
on the second orbital form an itinerant band»^ The two- 
orbital spin-only description would then amount to a ef- 
fective spin-localized description also of the itinerant elec- 
trons. Whether or not such a model can still correctly 
capture the low-energy magnetic properties is an inter- 
esting and important generic question. The aim of our 
work reported here is to establish some aspects of the 
phase diagram of the DSM more completely than in past 
studies,—"— as a starting point for investigating the low- 
energy physics in various phases and the nature of the 
quantum phase transitions between the different types of 
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FIG. 1: (Colored online) Six phases with coUinear order ap- 
pearing in the DSM. In the standard ferrormagnetic, Neel and 
stripe phases (upper row), the two spins on the same site are 
fully aligned. In the corresponding hidden-order phases (bot- 
tom row) there is no net on-site moment, due to the two spins 
on the same site being anti-parallel. 

ground states. 

Using a notation with a site index i e {!,..., A^} and 
an orbital index a £ {A, B\ on the S* = 1/2 spin opera- 
tors Sia , the DSM in the general form considered here is 
defined by the Hamiltonian 

^ = EE-^i^^s.^-s,,, (1) 

where and denote nearest- and next-nearest- 

neighbor pairs of sites, respectively, on the 2D square 
lattice. Here we count the A-B cross-terms only once, 
i.e., a > ^ in the sum over orbitals. This model has up 
until now been studied only in a very limited range of 
its parameter space.—"— We here report much more ex- 
tensive investigations, using a variational cluster-product 
wave-function ansatz. 

Our primary goal is to determine whether the DSM 
can provide a semi-quantitative, universal description of 
the magnetism of the pnictides and insights into their 
effective magnetic interactions. To address the issue of 
typically small but highly varying ordered moments in 
the (tt, 0) stripe phase, we take the view that the wide 
range of moments should originate from a sensitivity of 
the system to details of the effective spin couplings. It 
should therefore be possible to achieve small moments for 
a wide range of parameters in the DSM. The sensitivity 
should not require the proximity of a quantum-critical 
point, which seems to require too much fine-tuning to be 
generic. Note that the variations in the moments arc not 
only large among different members of the Fe pnictide 
family, but the moment of a given system can also be 
sensitive to pressure.— 

We also emphasize that the DSM is not just of interest 
in the context of the Fe pnictides, but also represents an 
intriguing and important quantum spin model in its own 



FIG. 2: (Color online) Illustration of average A (red), B 
(blue), and total spins (black) in a unit cell of the canted 
state. The canting angle 7 is obtained by rotating every sec- 
ond column (here the right one) with respect to the perfect 
Neel state (7 = 0), so that 7 = tt in the stripe state. 9 is the 
angle separating the A and B directions. 

right, with potential applications to other multi-orbital 
quantum magnets. It is interesting, e.g., to observe how 
the DSM evolves into an = 1 system with increasing 
Hund's coupling. More generally, we are interested in the 
quantum many-body states that can form in this system 
due to the two-orbital spin physics, beyond the known 
ground states of single-orbital S* = 1/2 systems and S = 
1 {Jh 00) systems. 

Using the variational cluster mean-field method with 
fully optimized cluster wave functions, as described in 
more detail below, in Sec. IIIH we find several ordered 
magnetic phases in the DSM, including ones previously 
identified and also some that had not been noted in the 
earlier studies of the double-spin model in a more limited 
parameter spacef^^"— To illustrate some of the complex- 
ity of the phase diagram, Fig. [3] shows simple collinear 
ordered and "hidden-order" phases that appear in the 
model in certain limits (that we will discuss in detail 
further below). Other, more interesting non-collinear 
phases with highly non-trivial quantum fluctuations ap- 
pear when tuning parameters between these limits. 

Being highly frustrated when at least some of the in- 
teractions are antiferromagnetic, the DSM is much more 
complicated and difficult to study reliably than the stan- 
dard Heisenberg model for the cuprates. Having two cou- 
pled spins per site, it is also much more complex than the 
5=1/2 J1-J2 Heisenberg model Spin- wave theoryi^ii^ 
and exact diagonalizationi^ have been used to study some 
regions of its parameter space. Retaining only Jh and 
the orbit-diagonal couplings jf^ and J^^, Manousakis 
discovered a canted phase with a 2 x 2 unit cell (illus- 
trated in Fig. [21)r^ in addition to the Neel and stripe 
states. He also suggested a non-magnetic phase for large 
Jh- Rodriguez and Rezayii^^i^ focused on completely 
different parameters, with only cross-orbital interactions 
J^^ and J^^ and weak Hund's coupling. States with 
"hidden order" then appear (Fig. [3]) where Neel or fer- 
romagnetic order is present but macroscopically invisible 
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because spins on the same site are anti-parallel (over- 
coming Jji) and lead to no net moment (if the effective 
(?- factors of the two electrons arc equal). These parame- 
ters do not seem appropriate for the pnictides, however, 
and the kinds of hidden order also seem inappropriate, 
because it does not correspond to the commonly observed 
(tt, 0) stripes in the Fc pnictides. 

The full phase diagram of the DSM ([T]) has not been 
studied so far. Even in the regions already studied, the 
methods used may not be completely reliable, due to 
the inherent difficulties with spin-wave theory when the 
quantum fluctuations are significant and the small lat- 
tices in exact diagonalization. It is therefore important 
to investigate the DSM for a wider range of parameters 
and using alternative methods. 



III. CLUSTER- VARIATIONAL METHOD 

The variational approach used here is based on fac- 
toring of the ground state of the system on the infinite 
lattice into a product of cluster states, 

i*)=ni^^)' 

c 

where '0c is the wave function of cluster c. We here use 
the standard basis of z-component spin eigenstates and 
write the state of a cluster of n sites as 

|^,)=^c<,|a), \a)^\Sl...S^^). (3) 

cr 

Imposing the constraint of all the cluster wave functions 
being identical corresponds to periodicity with a maximal 
unit cell given by the cluster used and mean-field clus- 
ter boundaries,— i.e., the interactions between different 
clusters factor into expectation value taken over two in- 
dependent clusters. Since all clusters have the same wave 
function, these expectation values just depend on a single 
cluster. In principle one may also extend the scheme to a 
product of different, collectively optimized cluster states 
(to allow for, e.g., magnetic order with larger periodic- 
ity) , but we will not consider this generalization here. 



A. Energy minimization 

Using clusters of size N = 2x2 and 2x4, we have fully 
optimized the wave function (which for the DSM has 4" 
coefficients) by minimizing the total energy, 

with respect to the parameters c^- To minimize E, we 
also compute its derivatives with respect to all the param- 
eters. It is useful to divide the energy into its diagonal 



(dia) and off-diagonal (off) parts, e.g., for the standard 
Heisenberg model (single-orbital) these are given by 

i?dia = Y.'^StS^^), (5) 

{i.j) 

Eos - i 5^(5+^7 +5.-5+), (6) 

and the generalization to the DSM energy is trivial. The 
general expressions for the diagonal and off-diagonal con- 
tributions to the energy derivatives are 

dEdia _ 2c„{Eaa — ij^dia) /-s 

dEoff _ J2t CtEctt - 'ic^EoS 

^ ~ — e:^ — ' 

where Ea-r = {'''\H\(t), with H is restricted to a single 
cluster and the mean-field decoupled inter-cluster inter- 
actions discussed above are included. 

The cluster-variational (CV) method outlined above 
is in practice very similar to the "hierarchical mean-field 
theory"—"— (producing identical results for systems such 
as the J1-J2 model^i), which in the past has been applied 
to several frustrated spin models. Our method of con- 
structing the wave function by direct energy minimiza- 
tion is different, however. We do not make any assump- 
tions regarding the form of the cluster wave function, 
treating all the coefficients in ([3]) as completely indepen- 
dent variational parameters. We normally do not imple- 
ment any symmetries (which would reduce the number of 
variational parameters) in order to not potentially miss 
any states. In principle symmetries can also be easily in- 
corporated in our scheme. Note that the magnetization 

n 

m. = Y.^SU + Sis) (9) 
1=1 

on the cluster is not conserved in general, but when 
studying the colincar states, such as the Ncel and stripe 
states, one can restric the wave function to the sector 
■niz = to reduce the number of parameters. In unre- 
stricted simulations the optimization does not necessarily 
find this subspace, because the spin-rotational symmetry 
is not broken explicitly (only spontaneously) and a state 
with m-y = in any direction 7 can be generated. All 
these solutions are degenerate and this is manifested in 
our calculations. 

We find the optimum using a combination of the 
steepest-decent method and a stochastic approach where 
only the signs of the derivatives are used,— with the lat- 
ter useful at the initial stage where it is important to 
avoid getting stuck in local minimums, and the latter 
working well at the final stage where the solution is close 
to optimal. For the models considered here, we have 
not encountered any significant problems with this op- 
timization scheme, even when the number of variational 
parameters is as large as 2^^ = 65536 (for a 4 x 4 single- 
orbital 5 = 1/2 Heisenberg model or a 4 x 2 cluster of 
the DSM). 
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B. Order parameters 

We here focus on the magnetic structure and define 
the Neel (Mi) and stripe (A/2) order parameters 
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where {xi,yi) are the integer-valued site coordinates. 
Naturally, these quantities as well can be expressed using 
only a single cluster in the CV scheme. 

The CV ansatz can break spin and lattice symmetries, 
but we do not impose how they are broken. We therefore 
compute the magnitudes of the full vector quantities ^TU\\ 
and In the case of stripes, AI2 corresponds to(7r,0) 
order (vertical stripes), but we also compute the (0, tt) or- 
der parameter and use whichever one (if any) is nonzero. 
In addition to the order parameters we also compute the 
canting angle 7 and the on-site separation angle 9 defined 
in Fig. [21 In calculations with 2x4 clusters we compute 
the magnetic properties on the central 2x2 plaquette. 



IV. RESULTS 

The magnetic structure within linear spin-wave the- 
ory is obtained by minimizing the classical energy, with 
quantum fluctuations only reducing (or possibly destroy- 
ing) that order. We will compare some of our results 
with spin wave theory. The CV method includes local 
quantum fluctuations within the clusters, and by using 
different cluster sizes we can check the stability of any 
states found. In practice we are of course limited to very 
small clusters, and cannot in general carry out completely 
unbiased extrapolations. The method nevertheless goes 
well beyond classical or standard mean-field approaches. 
The CV method should also lead to faster convergence 
with the cluster size in ordered phases (if this order is 
commensurate with the cluster structure) than standard 
exact diagonalization techniques. We here discuss sev- 
eral aspects of the phase diagram, starting with the case 
of Jh — > 00 (the spin-1 Heisenberg model) and then sev- 
eral other phases obtaining for finite Jh, including the 
parameter region we believe is appropriate for the pnic- 
tides. 



A. The S = 1 limit 



We first discuss the limit J// -> 00, where the DSM 
reduces to a spin-1 system. In the sub-space of parame- 
ters considered by Manousakis^^ we define g = J.^^/J^'^ 
(and all other parameters are zero). Classically, there is 
then a first-order transition between the Neel and stripe 
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FIG. 3: (Color online) (a) Neel (Mi) and stripe (M2) order 
parameters of the Jh = 00 DSM when the only non-zero 
parameter is g = J^^ I J\^ ■ The relevant order-parameter is 
Ml (Neel) for g > 0.5 and M2 (stripe) for g < 0.5. (b),(c) 
Dependence of the order parameters on g for three different 
Hund's couplings (using 2x2 clusters). Both Mi and M2 are 
nonzero in the canted phase. 



states at g ~ 1/2. As shown in Fig. EJa), within the 
CV approach with clusters of size 2x2 and 2x4, this 
transition moves only slightly, to g w 0.55. 

In linear spin-wave theory there is a narrow region at 
g « 0.5 where both ordered moments vanish (strictly 
speaking becoming negative, indicating a break-down of 
the approximation). Manousakis found that this behav- 
ior also persists for large but finite Hund's coupling, and, 
thus, that there is possibly a non-magnetic state inter- 
vening between the Neel and stripe phases. There is no 
explicit sign of such a state within the CV calculations, 
however. The Neel order drops significantly before the 
first-order stripe transition, and the drop increases with 
the cluster size. Such a drop is always expected in the 
presence of quantum fluctuations. In principle we cannot 
exclude that this behavior evolves with increasing cluster 
size into two independent transitions with an intervening 
non-magnetic state in between. 

In the = 1/2 J1-J2 model it is well known that a 
nonmagnetic state exists for 0.45 < 5 < 0.6J^"— This 
state is likely a valence-bond-solidr^ although some cal- 
culations instead have indicated a spin-liquid<^i^ A non- 
magnetic state in the DSM, if it exists, should then be 
the S = 1 version of such a state (which can exist^ for 
5 > 1/2 but is less studied than S = 1/2). Considering 
that the g-window of the nonmagnetic phase is already 
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FIG. 4: (Color online) Results for Manousakis' version of the 
DSM, where the only non-zero coupling is g — J2^/J^^ 
and Jh- (a) Evolution with g of the canting and separation 
angles (defined in Fig. [2]) at Jh = 10. (b) Boundaries in the 
{g, Jh) plane of the canted phase obtained in spin-wave theory 
(classical energy minimization) and within the CV method. 
In both (a) and (b) CV results are shown for both 2x2 and 
2x4 clusters. 



small in the 5 = 1/2 case, it should be very small indeed 
for S = 1 (since it must vanish in the classical, S" — > oo 
limit) and it is also possible that it is non-existent already 
for S = 1. The transition should then be first order on 
generic symmetry grounds. 

It is interesting to note that the location of the Neel 
to stripe transition is almost identical for the two clus- 
ter sizes, suggesting that the transition point is accu- 
rately converged. Moreover, the location also is close 
to the point obtained within a self-consistent harmonic 
approachi^ Note also that both ordered moments de- 
crease with increasing cluster size. Spin-wave theory is 
known to be accurate for 5 = (where the model reduces 
to the standard S* = 1 Heisenberg antiferromagnet).— The 
CV result is about 12% higher for these cluster sizes. 




FIG. 5; (Color online) The ground state energy of the DSM 
on the 2x2 cluster versus the coupling ratio g — J2^/Ji^ 
in Manousakis' versioni^ of the DSM at Jh = 10 in the = 
subspace (red) and the whole Hilbert space (black). The 
energies deviate from each other in the non-collinear canted 
state [with the canting angle shown for the same parameters 
in Ufa)] but coincide exactly in the coUinear states. 



7 and for JhIJ\^ = 10. The two different clusters 
give almost identical results. Note that the separation 
angle B is nonzero only inside the canted state, where 
< 7 < TT. The quantum phase transitions involving the 
canted state are clearly discontinuous for large Ju but 
become more smooth as Jh is reduced. 

The phase diagram in the (g, Jh) plane is shown in 
Fig. Hl^b) . The size of the canted phase is smaller than 
in spin-wave (classical) theory)^ but the fact that the 
boundaries are barely shifted when increasing the cluster 
size from 2 x 2 to 2 x 4 suggests that the state indeed exists 
and the computed phase boundaries should be stable. 

In a collinear state, such as those in Fig.[3l the ground 
state is in the sector of magnetization m,z = 0. The 
canted state is a non-collinear (co-planar), however, and, 
thus, its ground state should mix different magnetization 
sectors (regardless of the quantization axis chosen). We 
can confirm this explicitly by comparing the CV opti- 
mized energy in the full Hilbert space of the cluster with 
that obtained in calculations restricted to = 0. An 
example of these two energies as a function of g is shown 
in in Fig. [S] The region in which the two energies are 
different (with the one from the unrestricted calculations 
naturally being lower when they differ) coincides with 
the window in which the canting angle is not or tt. 



B. Canted state 

The canted state found within spin-wave theory by 
Manousakisi^ also exists in our approach for interme- 
diate g when Jjj is finite. Examples of the evolution of 
the Neel and stripe order parameters with g are shown 
in Fig. [3] for 2 X 2 clusters. Fig. SJa) shows the angles 



C. Modeling the Fe pnictides 

We now turn to our main objective of searching for 
parameters reproducing small and highly varying stripe- 
ordered moments of the Fe pnictides. The stripe moment 
tends to be large in Manousakis' version of the DSM (as 
seen in Fig.[3|) and one cannot expect the neglected model 
parameters to be very smalli^l We have investigated the 
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FIG. 6: (Color online) Stripe moments in the DSM with cou- 



plings Jh =2, Ji 



1, J2 



-0.7, and variable J2 



= 0.5, Jf " = 0.5, Jl'^ = 0.5, 
obtained with 2x2 clusters. 

-0.38. 



The separation angle 9 changes from tt to at J2 



full parameter space of the Hamiltonian ([1} extensively. 
In addition to the states discussed above and the hidden- 
order states studied in^^ we have also identified a hidden 
stripe state (see Fig. [3]) and a transition between it and 
the normal (tt, 0) stripes discussed above. 

Importantly, by considering couplings that are not 
symmetric with respect to the A and B spins, we can 
also achieve "partially hidden" stripe order, where the 
separation angle 9 = tt but the order parameter com- 
puted individually for the A and B sets of spins arc 
unequal. It is then possible to tune the net stripe mo- 
ment over a wide range of small to large values. This 
turns out to be the case only if the second cross-orbital 
coupling J^^ is ferromagnetic, in which case this cou- 
pling favors A and B stripes out-of-phase with respect to 
each other (and also further stabilize the stripes). The 
Hund's coupling must be small to moderate; up to 2 — 4 
times the largest exchange constant. This is perhaps 
smaller than expected in the Fe pnictideSfi^ but here it 
should be kept in mind that the bare Jh (in an elec- 
tronic model including the charge degrees of freedom) 
should be renormalized when reducing the itinerant de- 
scription to a localized-spin model. With localized spins, 
the effects of the Hund's coupling is naturally stronger, 
and, therefore, its bare value is renormalized down from 
its itinerant-electron value. All the other couplings can 
be antiferromagnetic, but a ferromagnetic J^^ also helps 
to stabilize stripes. 

Fig. [H] shows an example of the evolution of the sep- 
aration angle and the ordered moments as a function of 
J^^ when all other couplings are held constant at phys- 
ically reasonable^ values of order 1 (using J^^ as the 
energy unit). The canting angle jumps at ~ —0.38, 
corresponding to a transition between partially hidden 
stripes and normal stripes. There is no discontinuity in 
M2, however, because the individual moment of the A 
spins goes through zero at this point. The moment of 




FIG. 7: (Colored online) Phase diagram of the DS model in 
the subspace of varying J^^ and J^^ when all the other pa- 
rameters are held fixed at Jh ~ 5, Jj^^ = = 10 and 
Ji^ = J;^^ — 0. The three coUinear phases are ferromag- 
netic (F), stripe (S), and Neel (N). The hidden-order states 
are (where we have not drawn the boundaries that corre- 
spond to the coplanar partially hidden-ordered variants) hid- 
den (fully or partially) ferrormagnetic (HF), hidden (fully or 
partially) Neel, and hidden (fully or partially) stripes (HS). 
There are also three canted states, CI, C2, C3, with spin 
structure discussed in the main text. 



the B subsystem is almost constant at « 90% of its max- 
imal value 1/2 (and can be reduced by adjusting other 
couplings). At J^^ w —0.31 there is another interesting 
transition, where the stripe order jumps discontinuously 
with no qualitative change in the (tt, 0) structure. 



D. Other phases of the DSM 

The double-spin model studied here has six dimen- 
sionless parameters (coupling ratios) and a rich phase 
diagram in this space. To study phases beyond those 
discussed above, we have carried out systematic scans 
throughout the large parameter space. Here we show il- 
lustrative results obtained when some of the parameters 
are fixed within the subspace previously considered by 
Manousakis (and here in Sec. IIVB[) in the canted phase; 
Jh = 5, J^^ = J^^ = 10 and Jf ^ = J^^ = 0. We 
vary the cross-orbital couplings J^^ and to study 
how the canted state evolves. The resulting phase dia- 
gram is shown in Fig. [T] In addition to the canted state 



existing around Ji ~ J2 — (marked as CI in the 
figure) , all six coUinear phases illustrated in Fig. |3] are 
also present here. Moreover, there are two other kinds 
of canted phases here as well, marked in the figure as 
C2 and C3. These states interpolate between their adja- 
cent coUinear states, and correspond to the spin structure 
shown in Fig. [5] when switching the left and right spin of 
the bottom row (giving C2) or the top and bottom spin 
on the left column (giving C3). 
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V. SUMMARY AND DISCUSSION 

We have presented a comprehensive study of the two- 
orbital spin model introduced by Si and Abrahams Our 
results show that this system has a very rich phase dia- 
gram. In addition to the previously studied canted state 
interpolating between Neel and stripe antiferromagnetic 
states,— which we have shown to be stable against quan- 
tum fluctuations (although it is less extended than in 
spin- wave theory), there are also other canted states in- 
terpolating between all the three ordered collinear states 
we have studied (ferromagnetic, Neel, and stripe). More- 
over the previously studied hidden-order phases also have 
partially-hidden generalizations, in which the contribu- 
tions to the total ordered moment from the two orbitals 
do not cancel completely. The partially hidden stripe 
state is of particular interest in the context of Fe pnic- 
tides, reproducing the observed behavior of a sensitiv- 
ity of the magnitude of the moment to small variations 
in the parameters, without necessarily being close to a 
quantum-critical point. In this picture, the magnetism 
in the Fe pnictides is universally of the same type, but 
with different degree of cancellation of the moments re- 
siding on the two orbitals. 

One would expect antiferromagnetic couplings due to 
superexchange in the pnictides,— >2^i^ but the effective 
ferromagnetic cross-orbital couplings that we have found 
here in association with the partially hidden stripe state 
may appear due to the itineracy of the electrons. The 
superexchange was calculated in Ref. [13, but only as ef- 
fective Ji and J2 total-moment couplings, not separated 
into orbital-diagonal and cross contributions. A Hund's- 
like mechanism favoring parallel spins of A and B band 
electrons could come into play in the cross-orbital cou- 
plings. Such ferromagnetic couplings, generated as a sec- 
ondary effect by the on-site Hund's coupling, were in fact 



predicted very recently based on an electronic two-orbital 
model.— It would be interesting to also explore the role of 
effectively ferromagnetic couplings and partially hidden 
order also for the superconductivity of the pnictides<^ 

In principle it should be possible to detect the partially 
hidden aspect of the stripe order in NMR experiments, 
where different nuclei should probe the A and B spins 
differently and, thus, enable detection of their individual 
moments. This may require the hyperfine couplings to be 
determined more accurately than presently^ (including 
the different couplings to A and B spins). 

Beyond its relevance in the context of the pnictides, 
the DSM has a remarkably rich phase diagram, giving 
access to intriguing multi-orbital magnetic effects that 
arc not accessible with the standard single-orbital 5* = 
1/2 or 5 = 1 Heisenberg models. We have here pointed 
out several phases that had not been studied previously. 
It would also clearly be interested to carry out further 
studies to characterize these states and their excitations. 
There may still also exist other states in the large space 
of parameters of the DSM. 
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